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We theoretically study analytic-phase synchronization in strongly-competing oscillator systems. 
Using the example of composite-cavity modes coupled via a class-B laser active medium, we dis- 
cover that inherent chaotic phase synchronization can arise concurrently at two different chaotic 
£Sj ■ attractors, leading to bistable phase-synchronized solutions. In our example, the underlying mech- 

anism for bistability and inherent phase synchronization is population pulsation within the active 
medium. 
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Synchronization of interacting oscillators is encountered in many physical, chemical, and biological systems. Re- 
__, cently, considerable research has been devoted to understanding why and under what conditions chaotic synchro- 
^""^ ' nization is possible [1]. In chaotic synchronization, two (or more) chaotic oscillators adjust a given property of their 
r \ \ motion to a common behavior, as a result of coupling. Depending on the property under consideration, different types 
. ■ of chaotic synchronization were recognized; see [1] and references therein. Chaotic phase synchronization, which is 
Ph \ the subject of this paper, occurs when the phase difference of interacting chaotic oscillators remains bounded within 
i—i ■ the range of 2-7T for all time (phase locking) [2] . 

There are different approaches to achieving chaotic synchronization. A typical setup involves initially independent 
oscillators, where at least one is chaotic. Then, these oscillators (e.g. chaotic lasers) are made to interact with each 
fN| ' other, and the focus of existing studies is on either an unidirectional or a bidirectional coupling scheme [1-9] . Another 
^ , approach uses external modulation of a system of interacting oscillators to force chaotic synchronization, e.g. as in an 
experiment involving a modulated three- mode solid state laser [10]. Less understood and perhaps most interesting 
is inherent chaotic synchronization where the chaos and chaotic synchronization arise entirely from instabilities 
induced by the mutual coupling between the oscillators, without requiring chaos in uncoupled oscillators. Inherent 
chaotic synchronization was detected experimentally and studied numerically for lasers coupled at a distance [11-14]. 
However, the underlying physical mechanism remains unclear. Moreover, in all these studies, the chaotic-synchronized 
f^*) ' oscillators appear to involve a single chaotic attractor [1-14]. 

This paper investigates inherent chaotic phase synchronization within the scope of strongly-competing oscillators. 
Considering two composite-cavity modes coupled via class-B active medium [15] we are able to definitely identify 
(— | | the physical mechanism enabling phase synchronization to be an intrinsic behavior of this system. Furthermore, we 
discovered that the chaotic phase synchronization in strongly-competing oscillators exhibits bistability, a phenomenon 
that has not been reported previously. Also, we found that, in contrast to chaotic phase synchronization of oscillators 
that are independently chaotic, transition to inherent chaotic phase synchronization is not clear-cut. 

We consider a double-cavity laser, where cavity A of length L is coupled via a common mirror of transmission T to 
^ , cavity B of length L + dL. There are two approaches to model such a system [15]. (i) In the more phenomenological 
one, lasers are treated as individual oscillators and the coupling is introduced via ad hoc terms in the equations 
of motion for uncoupled lasers, (ii) In an alternate approach, called composite-cavity mode approach, the entire 
coupled-laser structure is treated as a single system. The lasing field is decomposed in terms of the eigenmodes 
extending over both cavities rather then the fields of individual cavities. In this picture, composite-cavity modes 
(CCMs), rather than individual lasers, are the interacting oscillators. The slowly varying dimensionless electric field 
amplitudes E n , and optical phases ip n associated with the n— th CCM evolve accordingly to [15,17] 

En = -jE n + C nnl x [CfcJ 1 + W A ) + c kn{^ + PN B )] cos(^ fc „) E h , (1) 

k 

i> n = Q n + C nnl x V [Ct (1 + pN A ) + C* (1 + [3N B )] sin^fen)!^, (2) 

k 

where 7 is the ratio between photon and population decay rates, /3 is the dimensionless gain coefficient, VL n is the 
dimensionless passive CCM frequency, and ipkn = V'fc — VVi [17]- In class-B lasers, the active-medium polarization 
decays much faster then the population and electric field. Then, the evolution of the dimensionless population N is 
governed by [17] 
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NA(B) = AA(B) - (N A{B) + 1) + - C™n ] (1 + f3N A{B) ) COs(Vnm) S m ^ n , (3) 

m,n 

where Aa(b) is the dimensionless excitation rate in cavity A(B). The T and cLL dependent Cnn 3 " 1 describes the 
overlap of the n-th CCM with the active gain medium in cavity A(B), and the nonlinearities due to optical coupling 
between the cavities [17]. We consider two CCMs, called symmetric (n = s) and antisymmetric (n = a) to refer to 
the relative optical phases between the fields in individual resonators [16]. For the calculations we use 7 = 2 x 10 11 
s-V(2 x 10 9 s- 1 ) = 100, (3 = 5.41, A = 2, and L = 280 [17]. 

The two sources of coupling (cross-saturation) between CCMs u n are spatial and spectral hole burning. In spatial 
hole burning, proportional to 

Cnm = G nn C mm = (4) 

i=A,B 

r r r L+dL r L+dL 
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competition arises because both CCMs deplete population at the same locations in the active medium [15]. In 
spectral hole burning, the electric field of one CCM saturates population at the frequency of the other CCM. 
One contribution to spectral hole burning comes from population pulsation, a result of nonlinear composite-mode 
interaction where the active-medium population acquires an oscillation at the intermode frequency [n ^ m terms in 
Eq. (3)]. This oscillation interacts with the n— th CCM to modify the active-medium polarization at the frequency of 
the m— th CCM. Consequently, more competition arises due to additional cross-saturation of the m— th CCM. Under 
appropriate conditions, the additional cross-saturation leads to strong competition resulting in bistability between 
stable stationary points [16,17]. 

Coupling between CCMs is adjusted with T and dL. For T = and dL ^ 0, C sa = and the two CCMs are 
uncoupled: each CCM becomes a mode of a different individual laser [16]. Recall, that we study inherent chaotic 
synchronization and the two CCMs are not chaotic at zero coupling. For |c?L| > 0, coupling increases with T [17]. 
As a result, the two interacting CCMs turn chaotic and exhibit transition to chaotic phase synchronization. 

Applying bifurcation continuation techniques [18] to Eqs. (1-3), we calculated saddle-node S, Hopf H, and torus T 
bifurcation curves in the parameter space (T, dL/X) [Fig. 1]. Curves S and H are tangent at saddle- node-Hopf points 
G where they change their type. Supercritical bifurcations of attractors are plotted as solid curves and subcritical 
bifurcations of unstable objects are plotted as dashed curves. Inside the lockband, which extends below the solid 
parts of S and H , composite modes are phase locked to operate at constant intensity and the same optical frequency. 
Moreover, there exist two stable stationary points in the {E a , E s , ip sa , Na, Nb} phase space. Depending on initial 
conditions, the system of coupled CCMs settles through strong competition to oscillate at an optical frequency that 
is 
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FIG. 2. (first column) Phase portraits and (second column) time evolution of the analytic-phase difference for the transition 
at dL/X = 0.55 x i<T 3 . From (a) to (f) T [1(T 2 ] =1.1, 0.97, 0.94, 0.93, 0.924, and 0.85. Refer to the dots a-f in Fig. 1. The 
time evolution of the analytic-phase difference is shown for the lower of the two coexisting attractors. 



near the frequency of cither symmetric or antisymmetric CCM. Each bifurcation curve in Fig. 1 denotes bifurcations 
of the two attractors [17]. 

Outside the lockband, optical-phase locking to a single optical frequency is lost and the modal intensities oscillate. 
Even when they are irregular, these oscillations may still exhibit certain types of synchronization when described with 
appropriate variables. Here, the appropriate quantity is the analytic phase ip k of a real signal x k (t) — E k (t) — (E k (t)), 

that is defined through x k (t) + ix k (t) = y 1 x\{t) + x\(t) exp[iip k (t)}, where (E%(t)) = lim T ^oo Jq [El(t)/T] dt, x k is 
the Hilbert transform of x k [19], and k — a,s. Analytic phase ip k is well defined if the trajectory in the complex plane 
{x k ,x k ) has one center of rotation. In the case of multiple centers of rotation, the signal needs to be decomposed 
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FIG. 3. (a) An example of a 2% jump in the analytic-phase difference [compare with Fig. 2 (d)] and (b) the corresponding 
chaotic oscillations of the coupled composite-cavity modes. 

intrinsic modes xu{t) — ^2jX^\t) such that the trajectory in the plane (x^\x^) has one center of rotation [20]. 

In this paper, we use ip^ whenever chaotic attractor has multiple centers of rotation. The analytic phase ifk of the 
modal intensity is distinctly different from the optical phase ipk- 

When the solid part of H is crossed from below to above in Fig. 1, locking of optical phases is lost. Each of the 
two coexisting stationary points turns unstable and gives rise to a stable periodic orbit. The relaxation oscillation, 
which is a characteristic of a class-B laser, become undamped. Now, both CCMs selfpulsate but their analytic phases 
are locked [Fig. 2 (a)]. Periodic oscillation can undergo further instabilities. Torus bifurcation curve T emerging 
from G marks the transition to quasiperiodic oscillation that involves two frequencies: the relaxation oscillation 
frequency and the CCM beatnote. Dynamics on the two tori also show analytic-phase synchronization [Fig. 2 (b)]. 
Decreasing T further leads to the break-up of coexisting tori into chaotic attractors. Despite chaotic oscillations in 
both CCMs, their analytic-phase difference remains within a range of 2ir for all time [Fig. 2 (c)]. Inherent chaotic 
phase synchronization is present simultaneously at two different chaotic attractors. Upon further decrease in the 
coupling, analytic-phase synchronization is lost [Fig. 2 (d)]. The analytic-phase difference for the lower (upper) 
chaotic attractor increases (decreases) in time exhibiting 27T jumps (Fig. 3). Next, the two chaotic attractors merge 
into one chaotic attractor [Fig. 2 (e)]. The analytic phases remain unlocked but the direction of 27r jumps alternates 
as a result of transitions between the two formerly bistable chaotic attractors. 

To see how the changes in analytic-phase synchronization of coupled chaotic CCMs come about, we plot in Fig. 4 
the maxima of E\ versus T. The simulation starts at higher T with two chaotic attractors, the lower and the upper. 
As T decreases, changes in the chaotic dynamics appear, including windows of periodicity called Arnold tongues. We 
find that analytic-phase synchronization is lost after transition through an Arnold tongue. At the right hand side of 
the widest Arnold tongue in Fig. 4(b), periodic orbit of high period is born and replaces the chaotic attractor. After 
some bifurcations, this orbit disappears to give place to a new chaotic attractor. However, the chaotic attractor that 
appears at the left hand side of this Arnold tongue [Fig. 5(a)] is significantly different from the chaotic attractor that 
disappears at the right hand side of this Arnold tongue [Fig. 5(b)]. The change(s) in the structure of the chaotic 
attractor (possibly due to homoclinic or heteroclinic tangency between stable and unstable manifolds of some saddle 
orbits) result in phase desynchronization. Analytic-phase difference for chaotic attractor at the left hand side of 
the Arnold tongue shows 2-7T phase jumps that become more frequent with decreasing T [Fig. 5(c)]. Similar changes 
take place in the upper chaotic attractor. Although the detailed bifurcation scenario inside the (narrow) Arnold 
tongue in Fig. 4(a) is different, the resulting effect is the same: phase-desynchronized chaotic attractor emerges. It is 
important to note that within desynchronized chaos we find Arnold tongues with dynamics that still show analytic- 
phase synchronization [arrow at the bottom of Fig 4(a)]. Therefore, transition to chaotic phase synchronization of 
independently stable oscillators is not as clear-cut as for oscillators that are independently chaotic [2]. 

As T decreases further, the lower chaotic attractor hits the basin boundary that separates the two chaotic attractors, 
and is destroyed [right arrow at the bottom of Fig 4(b)]. The trajectory settles to the upper chaotic attractor which 
then expands onto the remnants of the lower chaotic attractors [middle arrow at the bottom of Fig 4 (b)]. These 
two crises lead to a single chaotic attractor which is composed of the two formerly bistable attractors. Eventually, 
the chaotic attractor from Fig. 2 (e) hits the boundary that separates it from a periodic orbit [left arrow at the 
bottom of Fig 4(b)] and the trajectory settles to the periodic orbit with anti-phase dynamics [Fig. 2 (f)]. This orbit 
originates at the solid part of S and has a different origin than the periodic orbits in Fig. 2 (a). It arises due to 
unlocking of optical phases of the two CCMs [15,17]. Upon increasing T, this anti-phase oscillation coexists with 
in-phase oscillations, until it disappears in a saddle-node-of-periodic-orbit bifurcation at T <~ 1.1 x 10~ 2 . 

Different physical mechanisms are responsible for in-phase and out-of-phase synchronization. On the one hand, 
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FIG. 4. Bifurcations of the (a) upper and (b) lower chaotic attractor with decreasing T, leading to a loss of phase synchro- 
nization and subsequently to a single chaotic attractor. 




FIG. 5. (a-b) Change in the internal structure of the lower chaotic attractor due to transition through an Arnold tongue 
shown as a Poincare section defined by {Na = 0.025} [compare with Fig. 4(b)], and (c) the resulting loss of chaotic phase 
synchronization. 



when composite-mode detuning is too large for the optical phases to lock, but close to the relaxation oscillation 
frequency (T > 1.1 x 1CP 2 and outside lockband), population pulsation provides an inherent source of modulation 
that can excite self-sustained relaxation oscillation and force in-phasc dynamics of modal intensities. On the other 
hand, when CCM frequencies are close enough (T < 1.1 x 10~ 2 and outside lockband) the optical-phase locking terms 
become significant. Then, the system of two coupled CCMs alternates between the two states of nearly single-CCM 
operation, resulting in anti-phase dynamics. For intermediate conditions, in-phase and anti-phase dynamics can 
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coexist. Strong population pulsation effects are crucial for the inherent chaotic phase synchronization to occur. If 
they are neglected, neither bistability nor the instabilities leading to synchronization are present. 

In conclusion, we extended the usual analysis of chaotic phase synchronization to strongly-competing oscillators, 
where system dynamics are more diversified and less understood. In addition to previously reported results, the 
present analysis shows that chaotic phase synchronization can occur in coupled oscillators that are independently 
stable. This inherent chaotic phase synchronization results solely from the nonlinearities associated with strong 
mode competition within a saturable active medium. In the presence of population pulsation, strong competition 
causes the phenomenon to appear simultaneously at two different chaotic attractors, giving rise to bistable chaotic- 
phase-synchronized solutions. In contrast to phase synchronization of independently chaotic oscillators, transition to 
(inherent) phase synchronization of independently stable oscillators is not clear-cut: windows of phase-synchronized 
dynamics are found within phase-desynchronized chaos. 
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